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We explore the magnetic Feshbach spectra of ultracold ground-state Li-|-Er systems. Our calcula¬ 
tions predict many tunable resonances at fields below 1000 G that could be stably tuned in ultracold 
experiments. We show that Li+Er spectra are much less congested than those of systems involving 
heavier highly-magnetic atoms and exhibit non-chaotic properties. These features would facilitate 
identifying and addressing individual resonances. We derive a simple model for the mass-scaling 
shifting of low-field resonances that may simplify designing experiments with different Er bosonic 
isotopes. Our work establishes Li+Er as very promising systems for quantum simulation, precision 
measurements and the formation of polar paramagnetic molecules. 

PACS numbers: 34.50.Cx, 37.10.-x, 67.85.-d 


Ultracold species make it possible to build state-se¬ 
lected quantum systems with controllable interactions, 
which open the door to exploring fascinating phenom¬ 
ena. Among their many applications mm, ultracold sys¬ 
tems can be used as quantum simulators to study 

condensed-matter physics mm and quantum-controlled 
chemistry HUD], to develop quantum information de¬ 
vices mi HD and ultraprecise spectroscopy HMD]. 

Tunable Feshbach resonances |T6| are powerful tools 
to control the interaction and scattering properties of ul¬ 
tracold species, making many of these applications pos¬ 
sible. Moreover, they are essential in the most success¬ 
ful scheme to date to produce ultracold molecules: The 
magneto- or photoassociation of ultracold atoms mi¬ 
lls] followed by coherent |20fE5] transfer of the created 
molecules to their rovibrational ground states. Having 
the possibility to address and tune across selected Fesh¬ 
bach resonances is thus key in ultracold experiments. 

Recent advances in cooling highly-magnetic atoms 
such as Cr(^S) (SHISS], Dy(®I) jSDUlZ] andEr(^H) 


open exciting opportunities for tunability and control 
|30) . The interaction between these atoms, however, 
leads to highly congested Feshbach spectra with many 
overlapping resonances per gauss [3TH3D] . This makes it 
impractical to assign quantum labels to individual reso¬ 
nances and may be a challenge to interaction “tailoring” 
and molecule formation. 

In this paper, we study magnetic s-wave Feshbach res¬ 
onances in binary mixtures of ground-state Li atoms and 
bosonic Er isotopes. The Li+Er system may be spe¬ 
cially appealing for ultracold experiments in optical lat¬ 
tices: Dipolar species with tunable interactions are key 
to studying the effects of long-range anisotropies, quan- 
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turn magnetism, disorder and quantum collective behav- 
ior |ini l25l I34H36] . Very importantly, such Feshbach 
resonances may be used for magnetoassociation of LiEr 
molecules, starting from ground-state atoms in order to 
avoid limiting background losses m- Ground-state LiEr 
molecules have both magnetic and electric dipole mo¬ 
ments |38| . and may be controlled with applied electric 
and magnetic fields which further enhances their applica¬ 
bility H EH EDI |40] . In addition, Er is a heavy atom thus 
ultracold LiEr may be used to study the time variation of 
fundamental constants Ha mi, while the extreme mass 
imbalance in the system makes it specially well-suited for 
exploring Efimov physics [41] . 

We carried out coupled-channel calculations using the 
theory in Ref. [ID]- The Hamiltonian can be written 


n = 

2fj. dW 




+ Hu+HBr+U, ( 1 ) 


where pi is the reduced mass for the collision, R is the 
interatomic distance and L is the space-fixed operator 
for the end-over-end rotation. Hl, and R-Er describe the 
isolated atoms and are taken to be 


Ru = &F.Li*Li • SLi + {gsPBShi — ffiLiA^N^Li) ' B, 

R-Er = flEr^Er ' ^Er (SL^B^Er OSUBSEr) ' B. (2) 

Here, Ri and sli are the Li nuclear and electronic spin op¬ 
erators, while ^Er and sef denote the Er electronic orbital 
and spin operators (all bosonic Er isotopes have zero nu¬ 
clear spin); B is the external magnetic field, gs « 2, 
giEi and g/ 1 are the electron, Li nuclear and or¬ 
bital g factors, while /is and /xn are the Bohr and nu¬ 
clear magnetons. The hyperfine coupling constants for 
®’^Li (fep.Li) and the nuclear g factors were taken from 
Refs. I3DII33]- The spin-orbit coupling constant for Er, 
Ugj. = —1159.7215 X he cm“^, was calculated from the 
splitting between the two lowest Er states, ^Hg and ^Hg 
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|45| . assuming Russel-Saunders coupling. U describes all 
interactions between the atoms and includes the elec¬ 
tronic potential V and the direct dipolar interaction be¬ 
tween the atoms magnetic moments Hdip 111]- Following 
Krems et al. |46|, we decompose V into functions with 
well-defined total spin S and space-fixed spin projection 
Ms, which are then expanded in Legendre polynomials— 
assuming that Iet = 5 is conserved at all values of R. 
The expansion coefficients, V^{R) {k = 0,2,... ,21^^), 
are linear combinations of the Born-Oppenheimer poten¬ 
tials [1^ H5] : these can be split into isotropic and 
anisotropic terms, the latter depending only on the 
energy differences between Born-Oppenheimer states. 

The interaction between Li(^S) and Er(^H) gives rise 
to twelve electronic states: six states corresponding to 
|A| = 0,..., ^Er—the absolute value of the projection 
of the electronic orbital angular momentum onto the 
interatomic axis—for each S, S- = ^ and S+ = |. 
We calculated the short-range interaction energies for 
all Li(^S)+Er(^H) states using the complete active space 
self-consistent field method (CASSCF) implemented in 
MOLPRO m- The active space includes ls2s orbitals for 
Li and 4/6s6p orbitals for Er. We used the high-quality 
uncontracted aug-cc-pVQZ basis sets by Prascher et al. 
|48| for the Li atom. For Er, we used the quasirelativistic 
effective core potential by Dolg et al. jlHI (ECP28MWB) 
for the first 28 electrons, with uncontracted s and p shells, 
to which we added h functions with exponent 0.45. The 
basis was augmented with extra diffused functions using 
the even-tempered scheme in MOLPRO (47]. The Li+Er 
states lie relatively close in energy and special care is 
needed to avoid their mixing. We start our ab initio cal¬ 
culations at large R — bO qq, calculating the starting 
orbitals by merging those of the isolated atoms; this way, 
we obtain properly 22-fold degenerate states for Li+Er. 
In each following calculation, we take a step inwards in 
R and use the orbitals converged from the precedent ge¬ 
ometry as starting orbitals. We obtain potentials corre¬ 
sponding to pure A = 0,.. .±^Er states by controlling the 
A quantum number of the molecule. 

With this scheme, we get realistic short-range energy 
differences between all Li+Er potentials, and thus real¬ 
istic anisotropies I4 ^q, from our CASSCF calculations. 
The active space, however, is not sufficiently large for de¬ 
scribing the Li+Er dispersion accurately, which mainly 
affects the isotropic terms. This can be fixed by not¬ 
ing that the dispersion interaction in analogous systems 
is relatively spin-independent, with the spin-dependent 
exchange-dispersion energy being typically very small 
1^. We further assumed that the spin-independent dis¬ 
persion interaction in Li+Er is similar to that in Li+Yb, 
which is supported by the similarity of their isotropic 
dispersion coefficients (1594 a.u. for Li+Yb |2I] and 
1508 a.u. for Li+Er—see calculation details below). We 
used the Heisenberg spin-exchange model |52llH^ for the 
isotropic potentials, Vq{R) = Vo{R) — 2Jo{R)su-SEr, and 
replaced the spin-independent term Vg with the LiYb po¬ 
tential obtained by Zhang and coworkers using the best 
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FIG. 1. (Color online) Calculated interaction anisotropies in 
Li+Er for low (solid) and high (dotted) total electronic spin. 

S S 

ab initio methods available m- Jo = -[Po +-^o ]/2^+ 
depends on energy differences only and is thus recov¬ 
ered from our CASSCF calculations. We inter- and ex¬ 
trapolated all curves with the reproducing kernel 
Hilbert space (RKHS) approach of Ho and Rabitz |55| . 
Both isotropic potentials were constrained at long 
range to have Ce,o = 1 508 a.u., calculated from Tang’s 
combination rule with the values of the static po¬ 
larizability and dispersion coefficients for Li 2 m and 
Er 2 IHH]- The Vf anisotropies were constrained to have 
1^6,2 = 35.04 a.u., obtained from Cg.o and the anisotropic 
and isotropic polarizabilities of Er (c/.. Ref. [in]). We ne¬ 
glected the Van der Waals expansion in Vf and higher- 
order terms as they decay faster than 

Figure [l] shows the CASSCF interaction anisotropies 
Tfe^g for Li+Er. The magnitude of the spin-exchange in¬ 
teraction potential curves is qualitatively similar to re¬ 
cent calculations by Tomza on Li+Eu m, for which 
the spin-exchange interaction near equilibrium was about 
600 cm“^. The fact that spin-exchange in this system is 
much smaller than for alkali dimers is explained by sup¬ 
pression due to the outermost 6s^ shell of Er. The same 
mechanism reduces the anisotropies related to the elec¬ 
tronic orbital angular momentum: they are on the order 
of few hundreds of cm“^ near the Van der Waals min¬ 
imum for Vf terms and orders of magnitude lower for 
higher-order anisotropies. The mechanism of suppression 
of L-anisotropy was found earlier by Krems and cowork¬ 
ers for the He+transition metal systems m and for 
Yb+Tm by Buchachenko et al. (61]. This is in contrast 
to the Li+Yb(6s^6p^) system where the L-anisotropy is 
on the order of thousands of cm“^ m- 

We studied s-wave magnetic Feshbach resonances in 
the Li+Er systems using the molscat |62I[63] and field 
[51] packages. We used computational methods analo¬ 
gous to those in previous work on Li+Yb |37] and H+F 
[12] . The collision energy in our scattering calculations is 
fixed at 1 X fee nK. Convergence in the partial-wave ex¬ 
pansion was achieved with Lmax = 10. We consider cal¬ 
culations with both atoms in their lowest Zeeman state. 
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Magnetic field B (G) 

FIG. 2. (Color online) Magnetic field dependence of the scat¬ 
tering length a{B) for ground-state collisions of ^®®Er with: 
(a) ®Li, and (b) ^Li. 


Figure shows the predicted magnetic field depen¬ 
dence of the Li+^®®Er scattering length, depicting many 
Feshbach resonances of widths between 0.1 and 50 G 
that are very promising candidates for precise tuning and 
magnetoassociation. These resonances are immune to 
background losses since both species interact in their ab¬ 
solute ground state. We located 35/69 resonances below 
1 000 G for ®Li/^Li+^®®Er, and calculated mean densities 
of p6Li = 0.0364 G-i and = 0.0730 G-^: The com¬ 
bination of light and heavy species yields a much wider 
rovibrational spacing than for a heavy+heavy system, 
hence Li+Er Feshbach spectra are much less congested 
than others involving highly-magnetic atoms |31H33| . 

The couplings responsible for these resonances have 
been studied for analogous systems, c/.. Refs. EZllllllS]. 
These involve orbital- and/or spin-anisotropies from the 
electronic potentials, with the widest resonances due to 
states dominated by low L quantum numbers. The gen¬ 
eral trends in our calculations support this reasoning. 
The widest resonances in Fig.|^(AB > 10 G) correspond 
to states with over 55% and up to 45% compositions from 
L — 0 and 2, respectively, which are mainly coupled to 
the initial state by V/^q 2 terms. An intermediate res¬ 
onance “class” (1 < Ab < 5 G) arises from states with 
over 30% composition of T = 4, and less than 15% of 
L = 0, involving Vf coupling terms. Narrower reso¬ 
nances are due to states with non-negligible higher-order 
orbital excitations of the complex, L > 6 . Removing 
high-order anisotropies 14 ^ 4 , however, does not affect 
the overall resonance pattern since indirect mechanisms 
involving 14^4 terms also couple the resonant states to 
the continuum. The contribution from dipole-dipole in¬ 
teractions is very small compared to that of the potential 
terms and was found to be negligible, as for Li+Yb m- 

Studying the statistics of the calculated Feshbach spec¬ 
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FIG. 3. (Golor online) (a) Nearest-neighbor spacing distri¬ 
butions, and (b) number variances for the ®Li/^Li+^®®Er 
(solid red/blue) spectra. Poisson (dashed, green) and Wigner- 
Dyson (two-dashed, black) curves are added for comparison. 


tra is key to unraveling their most general and ro¬ 
bust properties. We performed an analysis similar to 
that of Frisch and coworkers [52], with added G-tests 
and Bayes’ logarithmic likelihood ratio (LLR) calcula¬ 
tions. Figure shows the nearest-neighbour spacing 
(NNS) distribution and the number variance for the cal¬ 
culated ®’^Li+^®®Er Feshbach spectra—predictions from 
the Poisson and Wigner-Dyson models are added for 
comparison. Table |T| summarizes all statistical tests per¬ 
formed and details on the fits to Brody NNS distribu¬ 
tions. All our tests indicate that Li+Er Feshbach spec¬ 
tra arise from weakly-interacting levels exhibiting very 
low spectra rigidity and level repulsion, with Poisson-like 
models providing a significantly better description of the 
spectra. These are all characteristics of non-chaotic spec¬ 
tra and make it possible to identify and tune selected 
individual Feshbach resonances for interaction tailoring 
and/or magnetoassociation. 

Erbium has five bosonic isotopes and the change in 
Li+Er reduced mass with respect to that of Li+^®®Er 
is about ±0.1% for ®’^Li. The small changes in /i, to¬ 
gether with the non-chaotic nature and relatively low 
densities of the spectra all support and would simplify 
predicting the Feshbach spectra for different Er isotopes. 
If SEy is the isotopic shift in the binding energy of a 
near-dissociation molecular state and (5/ires is the differ¬ 
ence between the magnetic moments of the molecule and 
free atoms (at resonance), the isotopic shift in resonance 
position, SB^es ~ SEy/Sfires, may be estimated for an 















































































































4 


TABLE I. Summary of statistical tests and/or model fitting to Li+ 


i66Er Feshbach spectra {x s-ncl G are reduced quantities). 


Lithium 

isotope 

Poisson 

’-test 

Wigner-Dyson 

(7-test 

Poisson Wigner-Dyson 

Bayes’ LLR 

log [P(a:|Hp)/P(a:|HwD)] 

Two-gap 

correlation 

Fit to Brody 

distribution 

®Li 

= 0.760 

X" = 2.53 

G = 0.982 

G = 2.08 

3.30 

rjr = 0.306 

p = 0.0738 


(p = 0.582) 

(p = 0.0435) 

(p = 0.436) 

(p = 0.0518) 



(MSE = 0.0101) 

"Li 

X^ = 1.40 

X^ > 10"° 

G = 0.760 

G = 8.06 

32.9 

Pr = 0.161 

p = 0.241 


(p = 0.246) 

(p < 10-®) 

(p = 0.654) 

(P< 10-“) 



(MSE = 0.000379) 


R ® potential to be |65l [66] 


cally representative for a much wider A range, except 


5B 


res 


r WKB _ ',2 

2-5Mres ^ '' 




where we approximated SEy « (d£'.u/d/r )|^^(^2 — 

Both the parameter Hqi « 3.4346/1^^1(700^^ and 
the WKB noninteger quantum number at dissociation, 
= (I/tt) arctan (1 — abgp/di) — 3/8—with di « 

0.47799 (2^i(76,o/^^) —are calculated for a reference 

system. As expected, Eq. (§ shows that Ji^res > 0 if 
1^2 > Ml) and viceversa (^Mres < 0 in our case). Eq. ([^ 
suggests that it may be possible to estimate the position 
of an equivalent resonance in a different isotopologue, 
while experiments may extract key information on back¬ 
ground scattering lengths, quantum numbers and non- 
adiabatic effects from isotopic resonance shifts. 

It is not yet possible to calculate quantitatively cor¬ 
rect ab-initio Li+Er potentials, hence the robustness of 
our conclusions need to be assessed for their possible de¬ 
pendence on these. The energy of the top vibrational 
state may be estimated from the Van der Waals coeffi¬ 
cient (76,0 s-nd reduced mass /r m to lie between 0 and 
13/i^/(/i7/^(3'w) ~ 6.3 X h GHz—where the characteris¬ 
tic length scale, i?vdw = ^(S/iC'e.o//!^)^^"^- The smallest 
difference between the magnetic moments of a scattering 
state and an L = 0-supported bound state is mb- The 
largest magnetic field at which resonances from such lev¬ 
els might occur is thus 6.3 x h GHz//rB ~ 5 000 G. Reso¬ 
nances resulting from larger magnetic moment differences 
may appear at much lower fields, however, 5 000 G is the 
smallest relevant range to study the dependence of our 
results on the interaction potentials. 

Figurefflshows a contour plot of the scattering length 
for ^Li+^Er as a function of the magnetic field, B, 
and a parameter A used to scale the isotropic poten¬ 
tials, Vq^ = AV)f. To keep a practical computational 
cost, the calculations neglect hyperfine terms and use 
Tmax = 6. The approximations preserve the main reso¬ 
nant features—compare with Fig.j^b)— while the broad¬ 
est resonances (A^ > 5 G) are most apparent at the 
grid’s resolution. Fig. reveals that there always ex¬ 
ist broad magnetic Feshbach resonances at relatively low 
fields. Additionally, it shows that A = 1 is statisti- 





Magnetic field B (G) 


FIG. 4. (Color online). Contour plot of the scattering length a 
for ^Li+^®®Er as a function of the magnetic field and a scaling 
parameter. Red/blue indicates the highest/lowest values. 


around 0.985 and 1.025 for which the background scatter¬ 
ing length becomes very large and enhances all resonance 
widths. 

This work provides robust theoretical evidence that 
low-field magnetic Feshbach resonances immune to back¬ 
ground losses exist for Li+Er, with widths Ab >0.1 well 
within current experimental resolution. Li+Er spectra 
are predicted to be non-chaotic, while remaining conve¬ 
niently dense, in contrast with other systems involving 
highly-magnetic atoms such as Er+Er [31]. The pre¬ 
dicted resonances may be resolved and independently 
addressed thus opening the door for precise tuning of 
Li+Er interactions and/or magnetoassociation into LiEr 
molecules. The characteristics of the spectra make it 
also possible to predict resonance positions for different 
isotopologues from measurements on a reference system, 
which would greatly simplify experiments with various Er 
isotopes and give key insight into non-adiabatic effects. 
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